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Abstract. Applying Baxter's method of the Q-operator to the set of Sekigu- 
chi's commuting partial differential operators we show that Jack polynomials 
( Xl ' ' ' ' ' Xn ) are e ig en f unc tions of a one-parameter family of integral op- 
erators Q z . The operators Q z are expressed in terms of the Dirichlet-Liouvillc 
• ra-dimensional beta integral. From a composition of n operators Q Zk we con- 

struct an integral operator S n factorising Jack polynomials into products of hy- 
pergeometric polynomials of one variable. The operator iS n admits a factorisa- 
■ tion described in terms of restricted Jack polynomials P^' ' 9 ' (xj, . . . , x^, 1, . . . , 

1). Using the operator Q z for 2 = we give a simple derivation of a previously 
known integral representation for Jack polynomials. 
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1. Introduction 

The method of Q-operator has been introduced by Rodney Baxter [3, 4] who used 
it to solve the eight- vertex model. Baxter's original Q-operator was a finite matrix, 
which commuted with the transfer-matrix of the model. Later on, Pasquier and 
Gaudin [31] discovered a connection between the classical Backhand transformation 
and a Q-operator for the periodic Toda lattice. Their Q-operator was a certain in- 
tegral operator with a simple kernel. Recent developments transformed the method 
into a unified approach used for solving quantum integrable systems, with several 
families having been extensively treated (see, e.g., [5, 23, 8, 19, 32, 34]). 

Let us take a quantum integrable system defined by n commuting linear partial 
differential (or difference) operators Hi in n variables, whose common eigenfunctions 
* A (x) =^x(xi,.-.,x n ), 

ff^ A (x) = ^(A)*a(x), (1.1) 
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form a basis in a Hilbert space 7i. The multi-index A = (Ai, . . . , A„) is a set of 
quantum numbers labelling spectrum hi{\) and eigcnfunctions %?\(x). 

By definition, a Q-operator Q z depends on a parameter z£C and satisfies two 
commutativity properties 

[Q Z1 ,Q Z2 ] = V Zl ,z 2 eC, (1.2) 
[Q Zl H,] = VzeC, Vi = l,...,n, (1.3) 
which imply that Q z can be diagonalised by the same basis functions ^a( x ) : 

[Qz*x](y) = nx(z)*x(y). (1.4) 

The third and the most important (as well as the characterising) property of 
a Q-operator is that its eigenvalues q\(z) satisfy a linear ordinary differential (or 
difference) equation with respect to the parameter z, 

w(z,^;{h i (X)})q x (z) = 0, (1.5) 



dz 

whose coefficients, apart from being functions in z, depend on the eigenvalues ftj(A) 
of the commuting operators Hi. Equation (1.5) is called Baxter's equation or, 
alternatively, the separation equation. 

An operator Q z for a quantum integrable system {Hi} is a quantization of a 
Backlund transformation for the corresponding classical integrable system. In the 
quasi-classical limit, its kernel turns into the generating function of the Backlund 
transformation, viewed as a canonical transform [31, 23]. 

A separating operator S n , by definition, factorises the basis functions v Pa(x) into 
functions of one variable: 

n 

S n : *A(Z1, ■ • -,X n ) >-> C\ (f)[ l ) (zt) (1.6) 

i=l 

c\ being some normalization constants. The partial functions <p x ^ (z{) are called 
separated functions and usually satisfy a differential or difference equation like (1.5). 
Notice that, in the quasi-classical limit, the operator S n turns into a separating 
canonical transformation for the corresponding Liouville integrable system. 

In [23], two of the authors noticed that any Q-operator gives rise to a family of 
separating operators for an integrable system. Indeed, given a Q-operator Q z and 
an arbitrary linear functional p on H, one constructs a product of n such operators 

Q. = Qz 1 ---Qz n (1-7) 

with the kernel Q z (y|x) as the convolution 

Q z (y|x)= J dti... J dt„_iQ Zl (y|ti)Q Za (ti|t a )...Qz n (t„_i|x), (1.8) 

and defines an operator S n p ^ = pQ z with the kernel 

<SW(z|x) = J dyp( y )Q z (y|x), (1.9) 

where p(y) is the generalised function corresponding to p. It then follows from (1.4) 
and (1.6)-(1.9) that S n P ^ is a family of separating operators (parameterised by the 
functional p) with 



[Zi)=q x (zi) and c A = J dy p(y)^x(y). (1.10) 

A clever choice of the 'projection function' p(y) can, in principle, simplify the 
structure of the integral operator S n p ^ ■ Our main motivation is to study suitable 
classes of Q-operators and functions p(y) for a variety of quantum integrable models 
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in order to construct simplest separating operators S n (1.6), which factorise special 
functions ^\{xi, . ■ . ,x n ). 

In the present paper, we apply the above ideas to the Calogcro-Suthcrland inte- 
grable model. In our case of study H is an appropriate Hilbert space closure of the 
space C[xi, . . . , x n ] Sn of symmetric polynomials in n variables, ^\(xi, . . . , Xn ) are 
Jack polynomials and A = (Ai, . . . , A„) is a partition. 

The problem of separation of variables for the Calogcro-Suthcrland model and 
for its q-analogue has been addressed by the authors in a series of publications 
since 1994. Many of the results presented below generalise those obtained from the 
studies of the few particle cases, n = 2,3,4, in [33, 20, 21, 18, 22, 24, 26, 27]. 

The paper is organised as follows. Section 2 fixes notation and collects the 
necessary background information on Jack polynomials. In Section 3 we study the 
properties of Dirichlct-Liouville multidimensional beta-integrals, which constitute 
our main analytical tool. All integral operators that we construct are expressed 
in terms of these integrals. In Section 4 we introduce a Q-operator Q z for the 
Calogero-Sutherland model and derive its properties. In particular, we prove that 
Q z is diagonal in the basis of Jack polynomials. Thus, the eigenvalue problem 
Q z ^fx = q\{z)^ \ can be considered as an integral equation for Jack polynomials. 
In Section 5 we obtain the eigenvalues q\(z) of Q z and give for them several explicit 
expressions, including the one in terms of the generalised hypergeometric series. 
In section 6 we construct the separating operator Sn^ and show that for a special 
choice of the functional p = p it splits into a product of n operators Ai, i = 1, . . . , n, 

5^)=^i ■■■A,, (l.H) 
such that each operator Ak acts only on k variables. Such a factorisation is an 
example of the factorised separation chain introduced in [17] for integrable models 
of sl(2) class. 

In Section 7 we describe another application of the Q-operator Q z . Setting 
z = we show that the operator Qq leads to an integral relation between Jack 
polynomials in n and n — 1 variables {xt}. Iterating this relation we produce an 
integral representation for Jack polynomials, which is equivalent to the integral 
representation due to Okounkov and Olshanski [28] who found it as a limit from a 
formula involving shifted Jack polynomials. 

2. Jack polynomials 

This section serves an introductory purpose as it fixes notation and collects the 
necessary background information on Jack polynomials. 

Jack polynomials P^ a \xi, . . . ,x n ) [12, 13, 14, 16, 35, 25] are orthogonal sym- 
metric polynomials, discovered by Henry Jack in 1970 as a one-parameter extension 
of the Schur functions s\{x\, . . . ,x n ). In the case a = 2, they reduce to zonal poly- 
nomials constructed by James in 1960 ([15]). We shall use the parameter g = 1/a, 
instead of a as in [25], i.e. use the notation P^ 9 \x) = p9~^ 9 \xi, . . . ,x n ). 

Let A = (Ai, A 2 , . . . , A„) e N™, Ai > A 2 > . . . > A„ > 0, be a partition of an 
arbitrary weight |A|, 

n 

|A|=5>. (2.1) 

i=i 

Hereafter, we shall use the weight notation | . | as in (2.1) for any finite sum of 
indexed variables. The dominance partial ordering ^ for two partitions \i and A is 
defined as follows: 

k k 

^ {|Ml = |A|; ^Mj<^A j; fc = l,...,n-l}. (2.2) 
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Usually (cf [25]), the length l(X) of a partition A is the number of non-zero parts 
A,. We shall consider the partitions A = (Ai, A2, . . . , A„) that are finite sequences 
of exact length n, including any zero parts. For instance, (6, 6, 2, 0) and (6, 0, 0, 0) 
are both partitions of the length 4, while (6,6,2) and (6) have lengths 3 and 1, 
respectively. The reason for this deviation from the standard definition is because 
we shall study Jack polynomials with the finite number of variables and charac- 
terise them as eigenfunctions of the complete set of n commuting partial differential 
operators {Hi} in n variables. In this picture, the parts Aj are the labels that pa- 
rameterise the spectrum and the eigenfunctions of the operators {Hi}. The number 
of parts, therefore, will always be equal to the number of the operators Hi (or to 
the number of the variables Xi). 

In the main text, for any partition A = (Ai, . . . , A„) of the length n, wc shall 
often use the notation Xij = Xi — Aj, in particular A^i+i = A^ — A,+i, i = 1, . . . , n, 
A n +i = 0. 

In the theory of symmetric polynomials, the number n of the variables Xi is often 
irrelevant because many results have been obtained for the case when this number 
is assumed to be large enough, for instance n > |A|, for the stabilised polynomials 
[35]. In this paper, the number n of the variables fixed finite number, the 

partition A is of the length n and has an arbitrary weight A|, which means that we 
shall mostly deal with the non-stabilised Jack polynomials. 

Let C(a, b, c, . . .) be the field of rational functions in indeterminates a, 6, c, . . . 
over the field C of complex numbers, <C(a, b,c, . . .)[x, y,z, . . .} be the ring of polyno- 
mials in the variables x, y, z, . . . with coefficients from C(a, b, c, . . .), and C[xi, . . . , 
x n ] Sn be the subring of symmetric polynomials. 

We need two standard bases in C[x] s ™ labelled by partitions A. For each a = 
(ai, . . . , a„) e N™ let us denote by x a the monomial 

— T Ql • • • T a ' 

The monomial symmetric functions m\ (x) are defined by 

ro A (x) = £V, (2.4) 

where the sum is taken over all distinct permutations v of A. The m^(x) form a 
basis in C[x] s ". 

For each r = 0, . . . , n the rth elementary symmetric function e r (x) is the sum of 
all products of r distinct variables Xi, so that eo(x) = 1 and 

e r (x) = ^ XilXi 2 ■ ■ -Xi r = TO(lr n-r)(x) (2.5) 

l<2i<...<2 T .<n 

for r = 1, ... ,n. For each partition A = (Ai, . . . , A„) define the polynomials E\(x) 
as 

£ A (x) - e t 1 - A2 (x) e ^- A3 (x)..-e r A l "(x). (2.6) 

The polynomials E\(x) also form a basis in C[x] s " and the transition matrix be- 
tween E\ and m\ is triangular with respect to the dominance partial ordering (see 
[25], Ch. 1, Sec. 2): 

^a(x) = ^2 °A/i"i M (x), (2.7) 

where e N and a\\ = 1. 

Jack polynomials P\ = P^ g \x) G C(5f)[x] 5 '™ are homogeneous of degree |A|: 

P| 1/9) (X!, . ..,x n ) = x^P^' 9 \ Xl /x n , . . . , X^/Xn, 1). (2.8) 



(2.3) 
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They form a basis in C(g)[x] s ™ . The transition matrix between the bases P\ and 
m\ is triangular with respect to the dominance partial ordering: 

P\ = v M m ^ v Xft {g) e C{g), v xx = 1. (2.9) 

The three bases, m\, E\ and Pa, are thus mutually triangular. 

The polynomials P^^(x) diagonalise the Sekiguchi partial differential operators 
[30] given by the generating function D(u; g): 



f_l)n(n-l)/2 

D(u;g) = A„(x) d6t 



D{u;g)P£' 9 \y) 



n 

Y[(\i + (n-i)g + u) 



.i=l 



P 



1/9) 



(X), 



(2.10) 
(2.11) 
(2.12) 



where 

A„(x) = det [x}" 1 ] " j=1 = J] (x, - x t ) 

l<i<j<n 

is the Vandcrmonde determinant. Expansion (2.9) and the fact that P\ are the 
eigcnfunctions of the Sekiguchi operators characterise Jack polynomials uniquely. 

Generating function D(u; g) is an operator polynomial in u whose coefficients 
give n commuting partial differential operators Hi, i = 1, . . . ,n: 



D( U ;s)= U " + f>"-%, Hl = ^-^g + ^i " 



i=l i=l 

Define the second order partial differential operator H as 



H 



n(n — l)(n — 2) 2 



g 2 + H 1 [H 1 -(n-l)g]-2H 2 



then 



H = Y( Xi — Y+gT^ ± ^ l U — -x-—). 
V 1 dxi J 4— ; Xi — Xj \ 1 dxi 3 dxi ) 



i<j 



It follows from (2.11) that 

n 

PP| 1/9) (x) = h(X)P { x 1/9 \ x ), h(X) = J>[Ai + ff (n+ 1 - 2i)]. 



(2.13) 

(2.14) 
(2.15) 

(2.16) 



Jack polynomials were applied in [6, 36, 29] to describe the excited states in the 
Calogero-Suthcrland model. Let us briefly remind some facts about this model. It 
describes a system of n quantum particles on a circle with coordinates < qi < n, 
i = 1, . . . , n. The quantum Hamiltonian and the operator of the total momentum 
are 



Hcs = - 



2 f-{ dqf 



E- 2 



9(9 - 1) 



sin {qi - qj) 



Pes = -i 



d 



3 = 1 J 



(2.17) 



The space of quantum states of the model is the complex Hilbert space L 2 (T n /S n ) 
of functions V&(q) on the torus T n = R n /nZ n 3 q, symmetric with respect to 
permutations of the coordinates g,. The scalar product in L 2 (T n /S n ) is defined as 
follows: 

(*,$) = !/ <%•••/ rfg„*(qT$(q). (2.18) 



The differential operators (2.17) are formally Hermitian with respect to the scalar 
product (2.18) when the parameter g is real. Most results of the present paper are 
established, however, for g > 0, which we shall assume hereafter. 
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The ground state of the model is given by the function 



w(q) = 



i<3 



(2.19) 



with the ground state energy E = ^g 2 (n 3 — n). The gauge transformation of the 
Hamiltonian Hcs (2-17) with the function w(q) gives 

w- 1 (q)oi/ cs o W (q) = ^H + E , (2.20) 

where the operator H is denned by (2.15) and the change of variables from q to 
x = e(q) is given by the map 

£ :(<?!,...,<?„) ~(e 2 ^,...,e 2i *"). (2.21) 

The eigenf unctions of the Hcs can therefore be chosen as the product of the ground 
state function and a Jack polynomial labelled by a partition A: 

* A (q)= W (q)pj 1/9) (e(q)). (2.22) 

The Hcrmitian form (2.18) corresponds to the following scalar product for any two 
symmetric polynomials /(x) and p (x) from C[x] 5 ™: 

(f,p)=[ dv Kv)] 2 /R^p( £ (v)), (2.23) 

where the integration domain, by symmetry, is 

T v = {v e R n /wZ n I Vj E (vj-uvj+i), j + n = j}. (2.24) 

The intervals (vj-\,Vj + i) above are understood as arcs on the circle R/ttZ such 
that Vj-i precedes Vj+\ with respect to the orientation on M. 

For g > 0, Jack polynomials p^ 1 ' 9 ^ (x) are orthogonal with respect to the scalar 
product (2.23). Note that p| 1/s) (x) = P ( x 1/9 \xT, ...,x^). 



3. DlRICHLET-LlOUVILLE INTEGRAL AND ITS MODIFICATIONS 

The Dirichlet integrals and their extension due to Liouville (see, e.g., [1]) are the 
most straightforward multidimensional generalisations of the beta-integral. They 
are used in computing volumes and surface areas of multidimensional ellipsoids. A 
special form of these integrals serves as the main building tool for our construction 
of the Q-operator for Jack polynomials. In what follows, the powers are functions 
of a positive argument £ > and are understood as = exp(aln£), with ln£ G K. 

Theorem 3.1 (Dirichlet-Liouville [1]). Define the domain C K" as 

^ = {€GR n |&>0,i = l,...,n}. (3.1) 

Then for SRai > 0, 

/ d* 5(\t\ - 1) ff 1 " 1 ' • • C"- 1 = r(a ^ (| ^ K) , (3-2) 

where a. — (ai, . . . , a n ), \a\ = a\ + ■ ■ ■ + a n , d£ = d£i ■ ■ • d£ n , and S is the Dirac 
delta- function. 

The proof can be found in [1]. We use two modifications of the integral (3.2) 
given in Theorems 3.2 and 3.3 below. 
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Theorem 3.2. Let y = (yi, . . . , y n ) be a fixed set of real positive parameters such 
that 

< Vi < V2 < ■ ■ ■ < Vn < oo. (3.3) 
Define the domain f2 x C M n by the inequalities 

fl x = {x e K n | < yi < x 1 < y 2 < . . . < x n _ x < y n < x n < oo}. (3.4) 

Let x = e„(x) = X\ ■ ■ ■ x n and y = e„(y) = y\ ■ ■ ■ y n . Then for jftc^ > and z > 1, 



f dx S(x zy) A„(x) f[& - 1 = (z- l)l°'l- 1 A n (y) 
J n* i=i 



r(|a|) 



(3.5) 



where A„ is i/ie Vandermonde determinant (2.12) and are defined as follows: 



(3.6) 



The other modification of the Dirichlct-Liouville integral is a trigonometric vari- 
ant of (3.5). 

Theorem 3.3. Given a fixed set of real parameters v = (v\, . . . , v n ) such that 

v x < v 2 < ... < v n < n + V!, (3.7) 
define the domain fi„ C 1" by the inequalities 

= i u € M™ | vi < m < v 2 < ■ . ■ < u n -i < v n < u n < tt + vx}. (3.8) 
Then for 7 e (0, n) and -Kc^ > 0, 



J 

JuT, 



du8{\u\ 



7) 



i<j 



where 



Li=l 

= (sin 7 )l a l- 1 \\-\nc, 

T = nLi sin K -vi) 



r(|a|) 



(3.9) 



(3.10) 



Both Theorems are proved by performing the respective changes of variables, 
from to {xi} (3.6) and from {^f} to {ui} (3.10), in the integral (3.2). 

Proof of Theorem 3.2. The first step is to rescale the variables & >— > — — 
in the integral (3.2) and to rewrite it in the form 



Vi(z-l) 



v*=l Vi ) i=l 



rriai 



(3.11) 



The next step is to prove that the relations (3.6) define a bijection of V£ onto 
f2 x . From the inequalities in (3.4) it follows immediately that tt x is mapped into 
V^. To prove the invertibility of this mapping consider the rational function 



i=i r Vi 

It is easy to see that rest= 3/i R(t) = — & and lim^oo i?(£) = 1, and therefore 

6 



fl(t) = i-£ 



(3.12) 



(3.13) 
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Given a set of positive £, > 0, the function R(t) is restored by the formula (3.13), 
and the variables Xi are then determined uniquely as its zeroes, taken in the in- 
creasing order. The standard continuity argument, using £j > 0, shows that all n 
zeroes Xj of R(t) are real and interspersed with the poles yj of R(t) and +oo or, 
equivalently, belong to the domain Q x . 

The bijectivity of the mapping £ i— > x having been thus established, its Jacobian 
is as follows: 



J = dct 



dxi 



det 



Xi yj 



A n (x) 

A„(y) 



>0, 



which was evaluated with the help of Cauchy's determinantal formula [7] 



det 



1 



1 — ciibj 



A„(a)A„(b) 



(3.14) 



(3.15) 



ij=i Ul^ii 1 -Oibj) ' 
Setting t = in (3.12) and (3.13) we get 1 + Ym=i €i/Vi = IIfc=i x k/yk, and hence, 



E 



- - (* - 1)) = (nw) * (n^ -*n^) ■ ( 3 - ie ) 

Using (3.14) and (3.16) we can now fulfil the change of variables £ i— > x in (3.11) 
and obtain (3.5). ■ 
The proof of Theorem 3.3 is similar to the above, but the structure of the map- 
ping £ T u (3.10) is, however, more involved compared to the rational case (3.6). 

Proposition 3.1. Define the domains , $J C K" as 

Vf = {£gR"| |£|<1, &>0;i = l,...,n} 

and 



= {u G nl I |u| < |v| + - }, $+ = {u G fi£ I |u| > |v| + - }. 
TTien i/ie change of variables u £ T maps the sets $J bijectively onto . 
Proof. Consider the function R(t) 

n 

R(t) = n 



(3.17) 
(3.18) 



sin(i — Ui) 
LL sm{t - Vi) 

Noting that reSi =Vi R(t) = —£[, R(t) can be rewritten in the form 

n 

R(t) = c -^([ cot(t - Vi). 
1=1 

Comparing the limits of the two expressions for R(t) at f — ► ±«oo one gets 

= sin(|u| - |v|) . 



(3.19) 



(3.20) 



(3.21) 



Suppose that u G fij, that is either u G $ u or u G $+. Then, from the 
inequalities (3.8) and the expressions (3.10) for , it follows that > 0. It also 
follows from (3.21) that £ T G (0, 1), for < |u| - |v| < tt. Thus, both domains $± 

are mapped into . Conversely, given £ T G , one obtains from the relation 
t 2 

+ £ = 1 (cf (3.21)) two possible values for the constant c in (3.20). The 

inequalities in (3.17) ensure that c 2 G (0, 1). Having chosen a solution for c, one 
reconstructs R(t) using the formula (3.20). The continuity argument, based on the 
positivity of £f, shows that the zeroes Ui of R(t) belong to the domain Oj. The 
positive c G (0, 1) correspond to u G <&~ and the negative c G (—1,0) to u G <£>„ . ■ 
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Proof of Theorem 3.3 starts with rescaling the variables i— ► -J^ in the 



integral (3.2) which transforms it into the formula 

eax-i . . . ea n -i _ ( sin7 )H-i r(ai)---r(q ra ) 



d€J(|€|-sin 7 ) ^ 



Sri 



r(|a|) 



(3.22) 



By the condition of the Theorem, 7 6 (0,7r), so sin7 S (0,1), and therefore the 
integration in (3.22) can be restricted to V£ . The Jacobian of the mapping £ T ^ u, 



J 1 



det 



det[£jcot(u 4 - Vj)]l j=1 , 



(3.23) 



can be calculated with the help of a trigonometric limit of the generalized (elliptic) 
Cauchy determinant (cf [11, 18]), 



det 



aim ~ Vj + C) 
a(xi - yj)a(() 



^(S + C) IL<j a ( x k - xi) a(y t - y k ) 



(3.24) 



ij=i CT (0 Ilfc,i ^ ~ w) 

where <r(a;) is the Weierstrass sigma-function and £ = ~ Di)- Take the 

limit of the infinite imaginary period, u>' — > 00, choose the real period w = ^ 
and put C = f j t nen cr ( x ) — * ex P (|s 2 ) sinx and (3.24) turns into the following 
trigonometric determinantal formula: 

det [oatfc - = cos( S ) ^ S ^ k ~ * l) ^ ~ ^ 



U k l sm(x k - yi) 



(3.25) 



Now we can evaluate the Jacobian J T in (3.23): 



J T = cos(|u|-|v|)n 



sin(uj — Ui) 
sin(vj — Vi) 



(3.26) 



Notice that J T is positive for u£$ u and negative for u e $+. 
Using (3.21), one transforms the (^-function in (3.22) as 

<J(|u|-|v|- 7 ) 



6(\£\ - sin 7) 



cos 7 



(3.27) 



Finally, combining (3.26) and (3.27) one fulfils the change of variables £ T 1— > u 
in (3.22) and obtains (3.9). When 7 G (0, n/2), the integration in (3.9) can be 
restricted to and when 7 <G (7r/2, it), the integration is over <!>+. Formula (3.9) 
is valid for any 7 £ (0, 1) by continuity. ■ 

4. Q-OPERATOR Qz 

In this section we construct a Q-operator Q z for the Sekiguchi partial differen- 
tial operators {Hi} (2.13). The operator Q z is described as an integral operator 
with a kernel expressed in terms of elementary functions and satisfies the defining 
conditions (1.2)— (1.4), proved in this section, and (1.5), proved in the next sec- 
tion. Thus, the eigenvalue problem Qz-Pj 1 ^ 9 ' = q\(z)P^ 9 ^ can be considered as 
an integral equation for Jack polynomials P^ 9 \ 

For g > 0, z > 1 and for an ordered set of n real positive parameters y = 

(yi, ■ ■ -,y n ), 

< yi < yi < ■ ■ ■ < y n < 00, (4.1) 

define the operator Q z , acting on symmetric polynomials p G C[x] s ™ , by the integral 



lQ*p](y) 



T(ng) (z-l) 1 -^ 



A n (y)f[»f- 1 



dx 5(x — zy) A„(x) 



nr 1 



p(x), (4.2) 
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using the notation of Theorem 3.2, notably 

_ IIfc=i( x fc ~ Vi) 

By virtue of the modification (3.5) of the Dirichlet-Liouville integral, Q z is nor- 
malised such that it sends p (x) = 1 to 1: 

Q»[l] = l. (4.4) 

The function [Q z p}(y) defined by the above integral for \/p G C[x] s ™ as a real- valued 
function of real arguments z and y is, in fact, a polynomial in z and a symmetric 
polynomial in y. 

Proposition 4.1. The operator Q z defined by (4.2) sends C[x] Sn into C(g)[y, z] s ™ . 

Proof. Let k — (fci, . . . , k n ) £ N n . Extend the domain of Q z from C[x] s " to 
C(y)[x]' s ", assuming that Q z acts trivially on the y- variables. In other words, wc 
can allow p in the formula (4.2) to depend on variables y in addition to x. In 
particular, we can set p = £ k . The integral (3.5) is well suited for calculating the 
action of Q z on the monomials £ k = ^ £ C(y)[x] s ". Substituting p = £ k 

into (4.2) and using the integral (3.5) for a>i = g + ki, we arrive at the following 
formula: 

/ 1 \ |k| n 

Q^w^^n^, (4.5) 

where 

(a) k = r( " + , fc) = a(a + 1) • ■ • (a + k - 1) Vfc e N, (4.6) 

is the Pochhammer symbol. 

The monomials | k form a basis in C[x] s ". To show this, it is sufficient to 
express the elementary symmetric functions ej(x), i = l,...,n (cf (2.5)) in terms 
of £, for any symmetric polynomial in x is uniquely expanded into the powers 
E\(yi) = e^~ A2 (x)e2 2 ~ A3 (x) • • -e A "(x) (see section 2). Introducing the generating 
function w x (t) for {ej(x)}, 

n n 

wx(t) = U(t xj) = t n + ^(-lyr-^x), (4.7) 

j=l i=l 

and using the formula (3.13) for the rational function R(t) — w x (t)/w y (t), we get 
w x (t) = w v {t)R(t) = w y (t) ■ (4-8) 

i=l 1 Vl 

Expanding (4.8) in powers of t, we obtain the expressions for {ei(x)} in terms of £: 

n 

ei (x) - e t (y) + (i [*-i(y)] Vj =o ■ ( 4 - 9 ) 
i=i 

We can give the following algorithm for calculating [Q z p}(y) for an arbitrary 
symmetric polynomial p e C[x] Sn : 

(1) express p(x) as a polynomial in {e^x)}, 

(2) substitute the expressions (4.9) for {ei(x)} in terms of {ej(y)} and 

(3) expand the resulting polynomial in monomials £ k , 

(4) replace each monomial £ k with the right hand side of (4.5). 
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The result [Q z p] (y) is obviously a polynomial in y and z and a rational function 
in g. Any permutation of the y-variables, a : y% i— > y ai , by (3.6), also sends & 
to . The right hand side of the (4.9) is therefore invariant under simultaneous 
permutations a : yi i— > y 0i and £j i— > Using (4.5), we conclude that the 

polynomial [Q z p](y) is symmetric in y. ■ 

Remark 4.1. Let ms mention another convenient variant of the above algorithm. 
After the substitution £j = r]jVj, the expression (4.9) for e^(x) becomes homogeneous 
in y: 

e i( x ) = ^2 ( 1 + r ?fci + ■ ■ ■ + Vk^Vki ■■■yk i - (4.10) 

l<fei<...<fci<n 

VKe can use (4.10) instead of (4.9) /or ffce step £ of the algorithm, and then use 
the substitution 

i i^lk ™ 

instead of (4.5), for the step 4- 

Remark 4.2. Notice also that the polynomiality of Q z £ k in z (4.5) ensures that 
the operator Q z , although initially defined for the real values z > 1, can now be 
analytically continued to the whole of the complex plane, z 6 C. 

Formula (4.2) is not the only way to realise Q z as an integral operator. An 
alternative description of Q z can be obtained from Theorem 3.3. 

For 7 e (0, 7r) and for an ordered set of real variables v = (vi, . . . , v n ) satisfying 

Vl < V 2 < ■ ■ ■ < V n < 7T + V\ , 

we define the operator Q 7 acting on a symmetric polynomial p £ C[x] s " by the 
integral 

KMM = r <"-"> <^>"" 9 



r(9)" TJsinfe-Di) 



i<j 



f du<y(|u|-|v|- 7 ) 



\ J s'm(u k - uj) 

j<k 



p(e(u)), (4.12) 



using the notation of Theorem 3.3 and the map e from (2.21). 
Proposition 4.2. Let 

x = e(u), y = £ (v), z = e('y). (4.13) 

Then [Q 7 p](y) G C(g)[y, z] s « and, moreover, Q 1 coincides with Q z as operators 
from C[x] 5 " to C[y] s ". 

Proof. Substituting p(e(u)) = (£ T ) into (4.12) and using the integral (3.9), 
one obtains the action 

which is similar to (4.5). Comparing expressions (3.6) for & and (3.10) for £^ and 
using substitutions (4.13), we obtain the relation 

0=ef-2M/^, (4.15) 
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which identifies the action (4.5) of Q z on the monomials £ k and the action (4.14) 

of Q 1 on the monomials (^ T ^j > an d therefore on polynomials p G C[x] s ™. Hence, 

Q z and Q 7 define two integral realisations of the same operator acting from C[x] s ™ 
toC[y] s ". ' ■ 

The integral formula (4.2) will be our main tool when working with the Q- 
operator. The trigonometric representation (4.12) will, however, be required for 
proving the conjugation properties of Q z (cf Proposition 4.4 below). In what follows, 
we shall identify the source space Cfx] 5 ™ and the target space C[y] s " of Q z and 
think of Q z as an operator in C[x] s ™. 

Proposition 4.3. Let e n be the multiplication operator acting in C[x] s ™ as e n ■ 
p(x) i— > e„(x)p(x). Then the following operator relation holds: 

Q z e n = ze n Q z . (4.16) 

Proof. For z > 1, the statement follows from the observation that Q z is realised 
as the integral operator (4.2) whose kernel contains S (e n (x) — ze n (y)). Then the 
result is extended by analyticity to z e C. ■ 

The following Theorem is the main result of this section, because it states the 
defining properties (1.2)-(1.4) of the introduced Q-operator. 

Theorem 4.1. Jack polynomials are the eigenf unctions of the operator Q z : 

lQ z Pi 1/9) ]( X ) = q x (z)P { x 1/a) ( X ), qx (z)eC(g)[z}. (4.17) 

The idea of the proof is to show that the matrix of Q z in the basis of Jack poly- 
nomials is both symmetric and triangular, and therefore diagonal. The symmetry 
is established by Proposition 4.4, the triangularity by Proposition 4.5. 

Proposition 4.4. The operator Q z possesses the conjugation property 

Ql = Qz (4.18) 

with respect to the Hermitian form (2.23). 

Proof. We shall use the representation of the Q-operator as Q 1 given by (4.12). 
To begin with, let us prove that 

(f,Q 1 p) = {Q*- 1 f,p) V/, P eC[x] s «, 7G(0,7r). (4.19) 

For our present purposes it is convenient to think that, similarly to (2.24), the 
variables u, v in (4.12) belong to the torus T n = M n /7rZ™ and all functions of u, v, 
including delta-function, are 7r-periodic. In particular, we can assume that all the 
sine functions in (4.12) take their arguments in (0, 7r) or, cquivalcntly, replace the 
sines by their absolute values. 

Keeping that in mind and using (2.23) and (4.12) we obtain 

(f,Q-yP)= [ dud V /(^v)yQ 7 ( V ,u)p( £ (u)), (4.20) 

where 

T UiV = {u,v e K"/7rZ n | uj G (vj,Vj+i), Vj G (uj-i, uj)} (4.21) 
(using the convention j + n = j) and 

Q 7 (v,u) = (sin 7 ) 1 -^^5(|u|-|v|- 7 ) 

x |sin(u fe - Vj) sin(ti fc - Uj)\ Y[ |sin(u fc - v^ 9 ^ 1 . (4.22) 

j<k j.k 
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Performing the following change of variables: 

Uj >—> u'j = Vj+i, vj >—> v'j = uj, (4.23) 
and observing that T u v — T u ', v > and Q 7 (v, u) = Q 7r _ 7 (v',u') we transform (4.20) 
into (Q w - 7 /,p). 

Having thus proved that Q 7 = Q w --y for 7 G (0,7r), we extend this result by 
analyticity to = Q n -^ for 7 e C. Due to the relation z = e 2 * 7 (cf (4.13)), the 
involution 7 7r — 7 corresponds to z 1— ► z, hence (4.18). ■ 

Proposition 4.5. T/ie matrix of the operator Q z between any two of the three 
bases, m\ (2.4), E\ (2.6) or P^ 9 \ is triangular with respect to the dominance 
partial ordering. For example, 

QzE x = ff)«V ( 4 - 24 ) 

Proof. Since the three bases, {m\, E\, P^ 9 ^}, are mutually triangular, it is 
sufficient to prove only the expansion (4.24). 

According to the algorithm explained in Remark 4.1, in order to calculate Q Z E\ 
we have to substitute into (2.6) the expressions (4.10) for ej in terms of r]j and yj, 
expand the brackets and transform the products of rjj by the formula (4.11). 

It follows from (4.10) that Q z maps E\ into a polynomial consisting of the same 
monomials in yj as E\ in Xj . The only difference is the presence of the variables 
r/j in the coefficients producing a nontrivial dependence on g. Since the expansion 
of E\ into m M is triangular (2.7), so is the expansion of Q Z E\. ■ 

Proof of Theorem 4.1. Let Q Xf , = (Q z P ( x 1/9] \P ( J~ /g) ). It follows from the 
triangularity property (Proposition 4.5) that ^ only for X h t 1 - On the 
other hand, from the conjugation formula (4.18) and the orthogonality of the basis 
p( 1/9} it follows that (Q Xfl ^ 0) => (Q m a ^ 0), therefore /i^A. However, A >r pt 
and \x >z A implies A = /x, hence (4.17). The eigenvalue q\(z) is a polynomial in z 
since [Q Z P^ /s) ](x) e C(g)[a;, z] 5 - by Proposition 4.1. ■ 

Remark 4.3. In the quasi- classical limit, g 00, the operator Q z turns into the 
corresponding Backlund transformation for the Calogero-Sutherland model ( see [23] 
for details). 

5. Eigenvalues of Q z 

Here we study the eigenvalues of the operator Q z , that is the polynomials q\(z) G 
C(g)[z] introduced in the previous section (cf (4.17)). 

First, we derive a representation for q\{z) as a finite multiple sum by making 
use of the explicit formula for the action of Q z on a polynomial basis. Then, we 
convert the obtained expression into another finite multiple sum, which results from 
a previously known formula (cf [20, 22]) for the separated polynomials q\(z) in terms 
of the hypergeometric function n F n -i{z). Finally, we write down the ordinary 
differential equation of order n for q\(z), which serves as the Baxter equation (1.5) 
for our Q-operator Q z , 

By Proposition 4.5, the eigenvalue q\(z) can be determined as the diagonal 
matrix element of Q z in the basis E\: 

Q z E x ^q x {z)E x + Y,bx^E^. (5.1) 

To calculate Q Z E\ we shall apply the algorithm described in Remark 4.1. Taking 
the diagonal element corresponds to replacing e^(x) by formula (4.10) and taking 
the coefficient at the monomial y A . Introducing the operator B g : C[»7] ^ C(g) 
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acting on the monomials by the formula (4.11) we arrive at the following ex- 
pression: 



q x {z) = B g 



Y[(l + m + ...+rj t ) x ^ 



(5.2) 



Theorem 5.1. The polynomial q\(z) £ C(.g)[z] is given by the multiple sum 

Ai,2 A n _i )Tl n—1 



fei=0 fe„_i=0 i=l 

Proof. Expanding the product in (5.2) we obtain 



in (-Ai,i + i)fc, (ig) kl +... + ki 

h\ ((i + l)g) kl+ ... +ki 



(5.3) 



Al,2 A n)Tl -|-l 

9aw=e- E 

fe 1= fc„=o 



n 



kj 



(z-l)l k l 
(«ff)|k| 



Ck 



(5.4) 



with 



Ck 



(ng) 



|k| 



B n 



(z-l)l k l ff 
(Di <8> • • • <g> £>„) 



3 = 1 



j=i 



(5.5) 



where the operators : C[^-] i— > C(g) are defined as follows : 77 - J > 
Lemma 5.1. 77ie coefficients Ck (5.5) are 

Ck = new + fci + . . . + k^ )k] = {n 9)M n ((3 -+i^x 



(5.6) 



The formula (5.3) is then obtained by substituting (5.6) into (5.4) and evaluating 
the sum over k n . ■ 

Proof of Lemma 5.1 is given by induction in n. When n = 1 formulae (5.5) 
and (5.6) give the same value {g)k x for c kl . It is sufficient to show that Ck lt ...,k n 
given by (5.5) satisfies the recurrence relation 



Ck u ...,k n = Cfci, • (ng + ki + . . . + k n -i) kn ■ 
Expanding the factor ((t]i + . . . + r/ n -i) + Vn) kn in (5.5) we get 

, n-1 

Cki,...,k. 



(5.7) 



(Di <g> ■ ■ ■ <g> £)„) 



l[(m + ■ ■ ■ + VjP )(vi + ■■■ + vn) h 
3=1 



E 



(£>i<g> •■•<8>£> n _i) 

\ /// / 

m=0 
ri-2 

3 = 1 

k 

c fci,...,/c T1 _2,^^-i+m {$)k n —m 



E 



m=0 



III 



Ck 



m=0 



III 



((n - l).g + fci + . . . + fc„-i) m (g)fc„-m- (5.6 
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To obtain (5.7) it remains to use the summation formula 

^2 ( k )( s )m(9)k-rn = (s+g)k (5.9) 
m=0 ^ ' 

for k = k n and s = (n — l)g + k\ + . . . + k n —i which is easily proved by comparing 
the coefficients with powers of t in the identity (1 — t)~ s (l — t)~ 9 = (1 — t)~ s ~ 9 . ■ 
In [20] (see also [22]) it was conjectured, and then proved for n = 2,3,4 (cf 
[20, 22, 26, 27]), that Jack polynomials P^ 9 \x\, . . . ,x n ) can be factorised into a 
product of the separated polynomials f\{z) defined in terms of the hypergeometric 
function n F n _\{z) as follows: 



f X (z) = Z X «(l-z) 1 - n <> n F n _! ( 



oi, o 2 ,o„ \ 
bi, . . .,b n -i J 

Oi = A„,i + 1 — (n-i + l)g, bi = cii + g. (5-11) 
Theorem 5.2. The function f\{z) is a polynomial given by the multiple sum 

Al,2 A n _l jn Ti—l . . 

fci=0 fe„_i=0 i=l 



fci! 



(l-(rt-j + l)g - Ai ]Tt ) fci+ .„ +fcn _ 1 

— 71 7 ^ \ — ^ ■ [p.iz) 

(1 - (n - ijg - A itn )k i +...+k n - 1 

Proof. Representation (5.12) for the separated polynomial f\(z) follows from 
the observation made in [22] that an explicit finite multiple sum for the polynomial 
given by (5.10) can be obtained by iterating Fox's formula [10], 



(7 \ 771 

ai,...,a p -i,b p -i+m \ _ sr^ ( , 
6i,...,6p_i ' Z J~h 



3=0 ■' 



p-1 / \ 
TT \ a i)] 

l\ Mi 



(ai+j, ...,a p -!+j 
p - 2 \b 1+ j,...,b p - 2 +j>' 



(5.13) 



which expresses the hypergeometric function p F p -i on the left, for any non-negative 
integer m, as a finite sum of lower series. 

Noticing that Oj+i = bi + Xi^+i for i = 1, . . . , n — 1 in (5.11) and using the 

binomial theorem, iF ( °; z ) = (1 — z)~ a , we apply (5.13) to (5.10) consecutively 



for p = n — 1, . . . , 2 and derive (5.12). ■ 
Theorem 5.3. T/ie polynomials f\{z) and q\{z) are proportional 

fx(z) = Pxqx(z), (5.14) 

where 

B, - TT /3 (l) /3 W - ((n ~' + 1)g)A »" f 5 15) 

Proof. We shall start with the expression (5.12) for f\{z). After the change of 
the variable z ■— > 1 — 2 we can rewrite (5.12) as 

/ A (l-z) = (l-z) A "^- 1) (z), (5.16) 
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where <p x \z; . . . , fc„_i) are defined recursively for i = 1, . . . , re — 1 by ip x °\z; 
ki, . . . , fc„_i) = 1 and 



^ ) (z;fc J+1 ,...,A : „_ 1 ) = z A -+ 1 ^ (— ) 

fe=o V z / 

(1 - (« - * + - Ai,„)fc i +...+k n _ 1 , (i-i). 



x v v- . „ ^-^ ( ^ }> (5 _ 17) 

(l-(n-t)ff-A i , n )fc i+ ... +fcn _ 1 



Lemma 5.2. 



T(g)T((n-i)g) 



<pX ) (z;k i+1 ,...,k n - 1 ) = (3 X 



rl 

X / d WU < ff- 1 (l- W )( n -09-l+A*+i.»-fci+i-----fc»-i( 1 _ w )AM + i $ (»-l)^ ) ( 51g ) 

Jo 

wftere $$^(-z) /or i = 0, . . . , n — 2 are defined recursively by & x °\z) = 1 and 

$f(z)= duu^^-u) 3 - 1 ^- zuf^+^^izu). (5.19) 
Jo 

Proof is done by induction. For i = 1, formula (5.17) gives 



^ (*; fa, . . . , fc„_ l} = (1 ^ g ' Al -" k+ - +fc 



n — 1 



(1 - (n - l).g - Xi,n)k 2 +...+k n - 1 



where 6*i = 1 — ng — Ai >ra + fc 2 + . . . + k n -\. Using the well-known identity for the 
terminating Gauss hypergeometric function, 

2Fl (-*> b ;z)=^(l- zY 2Fl ( 7; C - b ■ , (5.21) 



and then the integral representation (see, e.g. [9]), 

we arrive at the formula (5.18) for i = 1. 

Supposing that the statement is true for i — 1, we substitute the expression (5.18) 
for ^ in terms of ( f ) ^ 2 ' into (5.17) and get the integral 

i Px i) (z]k i+1 ,...,k n - 1 ) = z Xi - i+1 

rl 

x / d«w iff_1 (l -«)( n - i + 1 )»- 1 + A «.»- fc *+i— - fc »-i(i -zv)*'- 1 ''*!* -2 ^)^, 
Jo 

(5.23) 

where Si is the sum 

z-1 \ fel (-Ai, i+ i)fc 4 (1 - (n-i + l)g- \ i , n )k i +...+k ri - 1 



s i= 2 



fc s =0 



~' n V z ( 1 -' y )/ (! - ( n ~ i)9 ~ Ai, n )fe i +...+fc n _ 1 



(1 - (re - i + l)g - A i ,„) fc . +1+ ... +fcT ,_ 1 /-A^+i,^ (z - 1) . . . 
(l-(n-^-Ai, n ) fci+1+ ... +fen _ 1 2 X V ' 
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and 6i = 1 — (n — i + l)g — Aj jTJ + fcj+i + . . . + fc„_i . Using again the transformation 
(5.21) and then the integral representation (5.22), after some algebra we transform 
Si to the following form: 

(i) r((n-z + l)g) _ ^ 

^ r( 5 )r((n-z) 5 ) (1 ^ 

x / dtt 9 -\l - i)("-Ofl-i+A*+i.»-fc 4 +i-...-fc»-i [! _ y ( w + t (! _ „))] A -+\ 
Jo 

(5.25) 

Substituting Si from (5.25) into (5.23) we obtain an expression for ip^ as a double 
integral, in the variables v and t. Making the change of variables (v,t) i— > (v,u = 
v + t — vt) followed by (it, v) i-> = u/ti) we obtain for ip^ an expression as 

an integral in u and w. The integration domains are changed respectively from 
v, t G [0, 1] to < v < u < 1, and finally to u,w & [0, 1]. The result reads 



r(. 9 )r((n - 



^ l) (^;A; i+ i,...,fc n _i) = /3 A 



/ du« <fl_1 (l - u )("-*)»-i+Ai+i,»-**+i— -fc»-i(i _ zu ) A *..+i 
Jo 

/•l 

x / dww il - 1)9 - 1 (l - w) 9 - 1 ^ - zuw) Xi -^^ { l~ 2 \zuw). (5.26) 



The integral over w is evaluated with the help of (5.19) producing X \zu), and 
we recover formula (5.18), thus completing the induction step. ■ 
Iterating integrals (5.19) and using (5.16) we obtain the following representation 
for f\(l — z) as a multiple integral: 

h(l - z) = (1 - z) x -p x IM jf 1 dui . . . jf 1 

JJ wf-^l - ^^(l - *«< • • • «n-i) AM+1 . (5.27) 



n-1 

X 

i=l 



Expanding the binomials (1 — zui ■ ■ ■ u n _i) Ai i+1 we get a sum of products of one- 
dimensional beta-integrals in each Ui and, after evaluating them, we recover the 
expression (5.3) for q\(z), arriving thus at (5.14). ■ 

Remark 5.1. The eigenvalue q\(z) of the Q-operator Q z is thus proportional to 
the function f\(z) expressed by (5.10) as a hypergeometric polynomial, hence its 
Baxter (or separation) equation has the form (cf e.g., [9]) 

(4M B +t '- i )-^n(^ + »<)) [^< i -*>- i am] -o, 

\ i=l x ' i=l x ' / 

(5.28) 

where at = X n s + 1 — (n — i + l)g and hi = a, + g. 

6. Factorisation of Jack polynomials 

In this section, the integral operator Q z is used to construct a factorisation of Jack 
polynomials P^ 9 \x\, . . . , x n ) in n variables into a product of n polynomials in one 
variable. The general construction of the separating (factorising) operator S^ p \ 

S { n P) =PQ I1 ---Qz n , (6-1) 
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briefly described in the Introduction involves an arbitrary linear functional p : 
C[x] s " i ► C. Our present task is to make the integral operator si P ^ as simple as 
possible by choosing p = po in a special way. More specifically, we wish to factorise 
the integral operator Sn"^ into a string of operators Ak, 

S { n po) =PoQ Zl ---Qz n =Al---An, (6.2) 

possessing the property of dimension reduction in the sense that each operator Ak 
acts only on k variables. Any such factorisation will be called factorised separation 
chain or A- chain. 

Let us choose po to be the evaluation homomorphism that returns the value of 
a polynomial at the point x = (1, . . . , 1), 

PO :/->/*>/ = /(I,--., 1), (6.3) 
or, equivalently, the integral operator whose kernel is the product of delta-functions, 

p (x) = J] S(xj - 1), po.f = / dxpo(x)/(x) = /(l, . . . , 1). (6.4) 

3=1 J 

Having fixed the functional p — po, hereafter we shall omit p (or po) from the op- 
erator notation, writing S n to denote the corresponding separating operator. Two 
theorems below constitute the main result of this section and show how distin- 
guished is the chosen p (6.3). 

Theorem 6.1. The operator S n = PoQz x ■ ■ 'Qz n sends Jack polynomials into the 
product of separated polynomials q\(z) defined in Section 5, 



,(1/9) 

k=l 

the normalisation constant c\ being given by 



S n :P { x 1/9) (x)^c x l[q x (z k ), (6.5) 



l<i<j<n 



Proof. Apply the chain (6.5) of operators Q Zk to a Jack polynomial P^^(x) 
and use the fact (4.17) that P A Vs) (x) is an eigenfunction of Q Zk - Each factor Q Zk 
produces a desired factor q\(zk)- Finally, the evaluation homomorphism po (6.3) 
transforms Pj 1//9 '(x) into the normalisation constant c\ (6.6). For the evaluation 
of a Jack polynomial at x = (1, . . . , 1) see [35], [25, §6, (6.11')] or [28]. ■ 

Theorem 6.2. The operator S n : C[x] s ™ ^ C[z] s " ; corresponding to the evaluation 
homomorphism (6.3), can be decomposed into the product of n operators Ak, 

S n =A 1 ---A n , (6.7) 

A k :C[x 1 ,...,x k f k *->C[x 1 ,...,x k - 1 ] Sk - 1 ®C[z k ], k = l,...,n. (6.8) 

In particular, A\ : C[xi] >—> C[zi]. The operators Ak in (6.7) are assumed to act 
trivially (as unit operators) on the variables z. 

Furthermore, the operators Ak act on the restricted polynomials P^ 9 \x\, . . . , 
Xk, 1, . . . , 1) as follows: 

A k : Pi 1/9) (x!, . . . ,x k , 1, . . . , 1) h- P| 1/9) (^i, • ■ ■ ,a*_i, 1, . . . , 1) qx(z k ). (6.9) 

The factorisation (6.7)-(6.9) possesses the property of dimension reduction k 
k — 1. Such factorisations have been introduced in [17] when studying the inverse 
problem for integrable models of the sl(2) class and have been given the general 
name factorised separation chains. 
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Note that the formula (6.9) is consistent with Theorem 6.1. Indeed, applying the 
chain (6.7) of operators Ak to the polynomial P^ 9 \x\, . . . ,x n ) and using (6.9) 
we reproduce formula (6.5). The existence of a factorization of the form (6.7) is 
by no means trivial. One cannot simply define Ak by (6.9) because the restricted 
Jack polynomials P^ 9 \x\, . . . , Xk, 1, . . . , 1), as polynomials of k variables labelled 
by n-dimensional vector A, are not linearly independent. 

To prove Theorem 6.2 we shall describe a recursive procedure for constructing 
the operators Ak based on the defining formula S n — poQ Zl ■ • ■ Q Zn - The idea is to 
pull the homomorphism pa through the sequence of Q Zk 's transforming at each step 
Q Zk into Ak and producing a new homomorphism pk- Let idy^j for j < k denote 
the identity operator in C^ , . . . , Xk] Sk ~ i+1 ■ 

Proposition 6.1. There exist a unique sequence of operators Ak, 

A k : C^i,...,^]^ ' * C[x 1 ,...,x k -i} Sk - 1 ®C[z fc ], k = l,...,n, 
and a unique sequence of homomorphisms pk, 

pk ■ C[x k+ i, . . . ,x n ] Sn - k ^> C, k = 0, . . . ,n 

such that: 

(a) po : C[x] s " i— > C coincides with the evaluation (6.3), 

(b) the operators Ak are normalised by the condition Ak ■ 1 •— > 1, 

(c) the following relation holds for k = 0, . . . , n — 1; 

(id[i,fe] ®Pk)Qz k+l = A+i(id[i, fc+ i] ®Pk+i)- (6.10) 

The homomorphisms p k defined in this way act on polynomials f G C[xfc + i, . . . , 
x n ] Sn ~ k by restricting their n — k arguments to 1; 

Pk ■ f(xk+i, ■ ■ ■ ,x n ) ^ /(l, . . . ,1), fe = 0,...,n-l. (6.11) 

For k — n we have simply p n = 1. 

Proof of Theorem 6.2 follows then immediately. Applying recursively formula 

(6.10) to S n = poQ Zl Q Z2 ■ ■ ■ Q Zn we get S n = ^4i(id[ 1;1 ] ®pi)Q Z2 ■ ■ ■ Q Zn and so on, 

until we arrive at (6.7). Applying the identity (6.10) to a Jack polynomial P^ 9 ^ (x) 
and using (6.11) we get (6.9). ■ 

Proof of Proposition 6.1 is given by induction. For po the formula (6.11) holds 
by assumption. Assuming that for some k E [0, n — 1] the operators A\, . . . ,Ak 
and homomorphisms po, . ■ . ,Pk have been already constructed, let us consider the 
product (idhfci ®pk)Q Zk+1 and show that it factorises uniquely as .Afc+i(idfi ; fc_|_i] ® 
Pk+i)- Note that the operation (idfi^j ®pk) produces from a polynomial /(x) of n 
variables the polynomial f{x\, . . . , Xk, 1, . . . , 1) of k variables. As in Section 4, we 
will distinguish the variables y of the target space for the operator Q Zk+1 from the 
variables x of the source space. Given a polynomial / G C[x] s ", we have thus to 
take Q Zk+1 f G C[y] s ™ and to set yk+i = ■ ■ ■ = y n = 1- We can analyse the resulting 
expression in two ways based on two descriptions of Q Zk+1 given in Section 4: the 
one in terms of an integral operator and a pure algebraic one. 

Let us apply first the combinatorial prescription for evaluating Q Zk+1 on the 
elementary symmetric polynomials ej(x) given in Remark 4.1, formulae (4.10) and 

(4.11) . Recalling formula (4.7), we introduce the generating functions w x (t) = 
Ilj=i(* ~ x j) f° r e j( x ) an d' respectively, w y (t) = IIj=i(* — Vj) ^ or e j(y)- Consider 
now the relation (4.8), 

w x {t) =w y {t)-Y J ViVi ^r 1 , (6-12) 
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and set yk+i = ■ ■ ■ = Un = 1- The right hand side of (6.12) is then divisible by 
(t — l)™~ fe_1 , so is the left hand side, therefore n — k — 1 of the variables Xj are 
forced to take value 1. Due to the symmetry, the remaining free variables can be 
chosen as x\, . . . ,Xk+i- Note that from the above the statement of the Proposi- 
tion follows immediately: we have effectively transformed the operator expression 
(id[i,fc] ®Pk)Qz k+1 to the desired form ^4fe+i(id[i ; fc + i] <S>Pfc+i) with some operator 
Ak+\ acting on the variables X\, . . . , x^+i- 

We can go even further and find how Ak+i acts on polynomials in Xi, . . . , x k+ i. 
Let 

fc+i 

w x {t) = (t- l)™-*- 1 ^), w x (t) = H(t- x,) 

i=i 

and, respectively, 

fe 

w v (t) = (t-l) n - k - 1 w y (t), w y (t) = Y[(t- yj ), yo = l. 

j=o 

The relation (6.12) transforms then into 

w x {t) = w y (t) -^ViVi t 1 — . Va = Vk+i + ■ ■ ■ + Vn- (6.13) 

Expanding both sides of (6.13) in t we get the following prescription for calculating 
Ak+if- 

(1) express / £ C[x\, . . . , Xk+i] Sk+1 as a polynomial in {ej(xi, . . . , Xk+i)}, j = 
l,...,fc+l; 

(2) substitute for ej(x\, . . . , Xk+i) the expressions 

ej(xi, . . . ,Xk+i) = (1 + Vo + ■ ■ ■ + Vj-i)Vo ■ ■ • Vj-i + permutations (6.14) 

(the distinct permutations are taken simultaneously in yj and r/j); 

(3) expand the resulting polynomial in monomials ri™ ^™ 1 ■ ■ -ri™ k , m = (m , 
mi,...,TO fe ) eN fc+1 ; 

(4) replace each monomial in r]j according to the rule 

vo°vr ■ ■ ■ vr - (zfc+i ' i r°r' +T " fc((n ' fc)g)m ° nfow ( 6 - 15 ) 

(ng)m +...+m k jJi 

The last formula is obtained by expanding (i]k+i + - ■ ■+i]n) m " in the left hand side of 
(6.15), evaluating the result by formula (4.11) and using the following combinatorial 
identity for s = n — k, m = mo: 



U.(9)j a = {sg) m , (6.16) 



^ 7i 1 • • • 7 1 

ji + ...+j s =m J Js ' a = l 

which is easily proved by expanding the identity [(1 — t)~ 9 ] s = (1 — t)~ sg in powers 
of t. 

An alternative proof of Proposition 6.1 is based on the formula (4.2), presenting 
Q Zk+1 as an integral operator. Let us take the expression (4.2) for [Q Zk+1 .f](y) 
and restrict n — k of the arguments yj to the value 1. Since [Q Zk+1 f]{y) is a 
symmetric polynomial, it does not matter which of yj we choose to fix. Because of 
the inequalities in the definition of the domain £! x (3.4), 

< yi < xi < y 2 < ■ ■ ■ < y n < x n , 

the natural choice is to take the limit yi,t/2,- •• ,yn-k — > 2/0 = 1- The variables 
Xi, . . . ,x n -k-i are sandwiched between y's and, therefore, forced to tend to 1 as 
well. As yi, . . . ,2/n-fc tend to y = 1, the right hand side of (4.2) exhibits an 
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uncertainty 0/0. To resolve the uncertainty, let us set yj = 1 + evj , j = 1, . . . , n — k, 
and Xj = 1 + suj, j = 1, . . . ,n — k — 1, for some e > and Uj, Vj satisfying 

< vi < ui < . . . < u n - k -i < v n - k , 

and take the limit e — > 0. Let us also renumber the remaining arguments as 
Vn-k+j = Vj, j = 1, • • • , k, and x„- k +j-i = Xj, j = 1, . . . , k + 1. 

The components of the integrand of (4.2) are then transformed as follows: 



6 
-k+i 



I, 



i = 1, . . . , n 
i = 1, . . . k, 



where 



6 



A„(x) 



n"=i' ( u o - v J 

I\l=i k (Vj - Vi) 

Ujtlixj -m) 



1, . . . , n — k, 



0,...,k; 



rL:=° (yj - Vi) 1 

n+fc+l?n-fe-l A n _ fc _!(u) A k+1 (x) 



~ e 



1 1 n- 



A„_ fc (v) A fc+ i(y) 



A„(y) 

and 5 = (xi, . . . ,Xk+i), y = (yo, ■ ■ ■ ,yk)- 

Since dx\ ■ ■ ■ dx n ^ k -i ~ £ n ~ k ~ 1 dui ■ ■ ■ du n - k -\, the factors e cancel completely. 
Integration in the variables Uj produces a constant factor which is evaluated with 
the help of Proposition 7.2: 



i r 

A„_ fc (v) J V1 



dui 



dUn-k-l A„_fe_i(u) 



n — k 

nr 1 



r(s) 



n — k 



r((n - %) ' 

(6.17) 

The expression (id^^ ®Pk)Qz k+1 is presented in the desired form _4,fc+i(id[i.fc+i] 
(E>/9fe+i) with the operator ./U+i defined by the integral 



F(ng) 



(z k+1 - l) 1 ""* 



r((n - fc) 5 )rO/)* 



ife+i 



(v) 



m 

3 = 1 



/ 

JO, 



x / dx A k+1 (x) 5 (x-l ■ ■ ■ x k+1 - Zk+iyi ■ ■ ■ yk) Co 



'(n-fe)g-l 



nr 1 



/(SO. 

(6.18) 



where the integration domain f2 x is described by the inequalities 
< 1 = y < xi < yi < ■ ■ ■ < y k < x k +i < oo. 

Using the Dirichlet-Liouville integral (3.2), one can analyse the integral operator 
(6.18) in the same manner as the operator Q z has been analysed in section 4. In 
particular, one can prove an analogue of Proposition 4.1: A k +i sends symmetric 
polynomials into symmetric polynomials (with the dimension reduced by one) and 
its action on the elementary symmetric polynomials e 3 is described by formulae 
(6.14) and (6.15). ■ 

Notice that, for k = 0, the formula (6.18) produces A\ : f(x-\) i— > f(z\), whereas 
setting k = 1 in (6.9) we get 

Ai : P^ ,9 \x u 1, . . . , 1) ~ Pl 1/9) (1, . . . , 1) q x { Zl ) = c x q x ( Zl ). 
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Combining these two observations together, we arrive at the following remarkable 
expression for the separated polynomial q\(z) in terms of a restricted Jack polyno- 
mial: 

g A (z) = C - 1 F| 1/9) (z,l,...,l). (6.19) 

7. Integral representation 

In this section we shall apply our Q-operator to construct an integral represen- 
tation for Jack polynomials. In principle, such representation could be obtained 
by inverting the separating operator >S„, described in Theorem 6.2. However, the 
construction of S" 1 is a difficult and yet unsolved problem. Fortunately, there is 
another, more direct approach using the operator Q z at z = 0. 

Applying Proposition 4.3 to the case z = we conclude that Q nullifies the 
ideal e„(x)C[x] s ™ and therefore can be canonically defined on the factor space 
T = C[x\ Sn / e n (x)C[x] Sn . For the rest of this section we denote x = (x\, . . . , x n ), 
y = (yi, . . . ,y n ) and x' = (xi, . . . ,x n -i). Let V be the projection operator V : 
C[x] 5 " i ► CM 5 — 1 : p(xi, ...,x n ) i — ► p(xi, . . .,x n -i,0). Note that V : e;(x) ^ 
ej(x') for i = 1, . . . ,n — 1 and V : e„(x) 0. Since the products e^ 1 • • • f° rm 
a basis in T, the projector V provides a natural isomorphism T ~ Cfx'] 5 ™- 1 . We 
have thus come to the following conclusion. 

Proposition 7.1. There exists a unique operator Q' : C[x'] Sn - 1 i— > Cfy] 5 ™ such 
that 

Qo = Q' V. (7.1) 

Formula (7.1) provides a direct way to constructing a Jack polynomial in n vari- 
ables from a Jack polynomial in n — 1 variables. For any partition A = (Ai, . . . , A„) 
of length n define two partitions: A b = (Ai jn , A2, n , . . . , X n -i,n, 0) of length n and 
A 11 = (Ai jn , A2,„, • • • , A n _i,„) of length n — 1. Recall two important properties of 
Jack polynomials (see [25]): homogeneity, 

P\-/, 9 },\ n ( x l> ■ ■ ■ ' x n) — ( Xl ' ' ' x n) X ™ P\/^l..,\ n _ ln .o( x li ■ ■ ■ i x n), (7.2) 

and restriction to x n = 0, 

P^Lk-^A^ • • • ' ^-1'°) = < / „l,A n _ ll >i> ■ • ■ .^-i)- ( 7 - 3 ) 

Note that the equations (7.2) and (7.3) can be abbreviated as P^ g \-x) — [e„(x)] A " 

P^ /s) (x) and, respectively, V : PjJ /ff) (x) h-> p|J /s) (x'). 

By Theorem 4.1, the eigenvalue of Qo orL -Pa^'M ^ s gi yen by q\ (0). It fol- 
lows from Theorem 5.2 that, for small z, f\(z) = z Xn + 0(z Xn+1 ) and then from 
Theorem 5.3 that 

q x (z)=p- 1 z x "(l + 0(z)), z^0, (7.4) 

where (3\ is defined by (5.15). In particular, q\t>(0) = /J^ 1 . 

Applying the operator equality (7.1) to the polynomial P^ 9 \x) and using (7.3) 

and q A b(0) = we obtain fc 1 P§ /9) (*) = [g o P{J /s) ](x). Finally, we use (7.2) 
and arrive at the formula 

if /9) (x) = (3 X [e„(x)] A " [g p(J /9) ](x), (7.5) 

expressing a Jack polynomial in n variables in terms of a Jack polynomial in n — 1 
variables. 

To give the formula (7.5) more flesh, we need to realise Qo as an integral operator. 
However, formula (4.2), which we used to define Q z , is valid only for z > 1, so that 
the case z = requires a special consideration. 
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Proposition 7.2. Given a set of real parameters y = (t/i, . . . , y n ) satisfying in- 
equalities (3.3), define the domain C W 1 ^ 1 by the inequalities 

Q' x = {x' e K™- 1 I < yi < x x < y 2 < . . . < a;„_i < y n < oo}. 

Then for Sftcq > 

, A„_!(X') 



(7.6) 



/ 



dx! 



A„(y) 



.»=i 



r(ai) ■ ■ ■T(a n ) 
r(|a|) 



where & > are defined as 



11™=! (j/fc - ' 



i = 1,. .. ,n, 



E6 = 1 - 



(7.7) 



(7.8) 



i=l 



The proof follows closely that of Theorem 3.2. Introducing functions w x (t) — 
riiLi — x «) an( i w y(t) — n"=i(^ — 2/0 we construct the generating function 
R(t) — w x (t)/w y (t) = Y^i=i &/(* — Vi) f° r * ne cnan g e of variables (x\, . . . ,x n -i) 
(Ci) ■ ■ • ttn-i)- The argument, based on studying R(t), shows that the domain Sl' x 
is mapped bijectively onto V£, 

V? = {( , 6R n - 1 |&>0, i = l,...,n-l, 6 + ... + |„_i<l}. (7.9) 
The corresponding Jacobian is given by det [9^-/9a;fc]™ . =1 = A„_i(x')/A„(y). 

Making the change of variables x' £' in the integral (7.7) we come exactly to 
the Dirichlet-Liouvillc integral (3.2) with the only difference that the integration 

variable is eliminated by resolving explicitly the constraint f =1. 
Theorem 7.1. The integral operator Q defined by 



[Qp](y) 



r(»g) 



A n -i(xQ 
A„(y) 



nr 1 



.4=1 



p(x') 



(7.10) 



sends Clx.'} 8 '™- 1 into C[y] 5 ™ and coincides with the operator Q' defined by (7.1). 

Proof. Repeating almost word by word the proof of Proposition 4.1 and using 
formula (7.7), we find the analogue of the formula (4.5) Vk = (fci, . . . , k n ) E N™, 

(s)fci • • ' (ff)k„ 



(7.11) 



("ff)|k| 

The action of Q on the polynomials ej(x') is obtained then from the analogue of 
the formula (4.9), 



e *( x ') = l e i-i(y)]i/ 3 =o' 



(7.12) 



i=i 



which is obtained, in turn, by expanding the equality w x (t) — J27=i ii w y{t)/{t~yi) 
in powers of t. 

It now remains to identify the action of Q and Q' on Cfx'] 5 " -1 . Setting x„ = 
in the expressions (3.6) for & and comparing the result with the expressions (7.8), 
we obtain 

Vii = -ytii. (7.13) 
Setting 2; = in the right hand side of (4.5) and using (7.13) we obtain exactly the 
right hand side of (7.11), which proves that Q a = QV. ■ 
We can regard the formula (7.5) as an integral relation, which reduces a Jack 
polynomial in n variables to a Jack polynomial in n — 1 variables 1 . Starting with 



1 or, vice versa, builds up the number of variables, starting from the unit function 
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P^ 9 ^ (x\ ,x n ) and having iterated the formula (7.5) n — 1 times, we go down to 

the P^J 9 \x\) — x^ 1 and, as a result, get an explicit integral representation for a 
Jack polynomial in n variables. To describe it explicitly, we introduce a triangular 
matrix Xij, 1 < i < j < n of variables. At the first step, the variables 
i = 1, . . . , n, are the external ones and x^„_i, i = 1, . . . ,n — 1, are the integration 
variables. Then the equation (7.5) takes the following form: 



P (1/s Vri r 1 



n 




W x i,n 









n 



((n-i + l)g) Ai „ 
L = l (( n - »)5)a*,„ 

dx n —\ iT i—\ 



1 • • ' 



Yli<j{ x j,n—1 Xi^ n -i) 
Il«<7 — x i,n) 



' &n — 1 ,« 
n n— 1 



i=i fc=i 



9-1 



X P, 



\ n _ 1 n (^l,n— 1) • • • j l,n— l)* 



(7.14) 



Iterating this formula, we obtain the integral representation for Jack polynomials 



n 



n v n 

j=n— 1,1 i=l,j 



i<j 

dxi 

Xi,3 + 1 



A 3 rf ((j + 1-%)A., 

rife=l( a; ij — x k,j-lj 



nn 



n-l J 



nn 

i=l i=l 



9-1 



(7.15) 



It is implied in (7.15) that the integrals J dxij are ordered corresponding to the 
order of indices in the products. 

As mentioned before, formula (7.14) first appeared in [28], where its proof was 
based on the fact that Jack polynomials can be obtained as a limit from the shifted 
Jack polynomials. One can regard the integral equation (4.17) in section 4 as the 
one-parameter generalisation of (7.14). 



8. Concluding remarks 

Our construction of the Q-operator relies on the properties of Jack polynomials. 
Nevertheless, several of the obtained results are of importance for the general the- 
ory of quantum integrable models. It concerns, first of all, the factorised separation 
chain, constructed in Section 6. Similar factorisations of separating operators were 
already observed for other integrable models [17] and, apparently, they are mani- 
festations of some general pattern, which is yet to be fully understood. 

Notice also the intriguing role of the restricted Jack polynomials revealed in 
formulae like (6.9) or (6.19). It would be interesting to find analogous formulae for 
other integrable models. The properties of restricted Jack polynomials themselves 
deserve a further study. 
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Finally, we should mention another integral representation known for Jack poly- 
nomials [2], where, in contrast with (7.15), the number of integrations grows with 
the weight of A. 
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